Gate-controlled spin polarized current in ferromagnetic single electron transistors 
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We theoretically investigate spin dependent trans- 

port in ferromagnetic/normal metal/ferromagnetic single electron transistors by applying master 
equation calculations using a two dimensional space of states involving spin and charge degrees of 
freedom. When the magnetizations of ferromagnetic leads are in anti-parallel alignment, the spins 
accumulate in the island and a difference of chemical potentials of the two spins is built up. This 
shift in chemical potential acts as charge offset in the island and alternates the gate dependence of 
spin current. Taking advantage of this effect, one can control the polarization of current up to the 
04 . polarization of lead by tuning gate voltages. 



PACS numbers: 72.25.-b, 73.23.Hk, 85.75.-d, 85.35.Gv 

Ferromagnetic single electron transistor (SET) has been an interesting system shown-to exhibit novel phenom- 
ena with an interplay between spin and charge. Recently, Chen et a» and Ono et aa succeeded in fabricating 



small double junctions containing metal or superconductors sandwiched between two ferromagnets. In their exper- 
iments, enhanced tunneling magnetoresistance(TMR), magneto Coulomb oscillations and spin accumulations were 
observed. On the other hand, theories of ferromagnet/ferromagnet/ferromagnet(F/F/F) and ferromagnet/normal 
metal/ferromagnet(F/N/F) SET's based on transition rate and master equation formalism were developed to de- 
rive bias- voltage and gate- voltage dependent TMR in both sequential and strong tunneling regimes The pio- 
neer experiment conducted by Johnson and Silsbee demonstrated the importance of spin accumulation effect on 
ferromagnet-normal metal systemsij For F/N/F double junctions, the spin accumulation is predicted to occur when 
two ferromagHgjUleads are in anti-parallel alignment, which would lead to a new origin of TMR in contrast to that of 
F/F/F cascs.oQu In this study, we investigate the spin accumulation and related phenomena in F/N/F SET under 
O \ the influence of gate charge. 

The spin dependent transport in a ferromagnet is usually described, by the relative difference of the majority and 
—i | minority spins of conduction electrons, denoted as polarization P.E3 Under the condition that spins do not flip, 
^ i the transport current in F/N/F double tunnel junctions can be separated into two channels labeled as upspin and 
downspin which, throughout this article, are assumed to be contributed by majority and minority spins, respectively, 
in the source ferromagnet. For example, when lead magnetizations are in anti-parallel alignment, the upspin channel 
has a larger tunneling rate for the source junction than for the drain junction. In this case and in the steady 
condition, the upspin chemical potential in central electrode rises to balance the spin's incoming and outgoing rates, 
and the chemical potential of the downspin would decrease by the same amount. This shift in spin chemical potential, 
denoted as A/x-|-m, is predicted to be ±PeV&/2, in which Vt> is the applied bias voltage (see inset of Fig. 1) and P is 
the polarization of the two leads. Therefore, the net spin in normal metal becomes nonzero and this effect is known 
C$ ' as the spin accumulation (or spin imbalance). 

A gate voltage V g is applied to turn on and off the charge transport in SET by tuning the electrostatic potential of 
the island. When SET is symmetrically biased with V g =0, an energy cost of roughly the charging energy Eq is required 
for adding or removing an excess charge in the island, and the charge transport is blockaded provided that temperature 
is low, ksT -C Ec- When the gate voltage is tuned with two adjacent charge states becoming energetically degenerate, 
electrons can enter or leave the island without extra energy cost, producing sequential tunneling current. Ideally, the 
current is at a minimum and maximum, respectively, for V g = and e/2C g and can be modulated periodically with 
. £h ' a period of AV g — e/C g ; here C g is the island-to-gate capacitance. These two gate voltages are thus referred to as 
minimal and maximal gate voltages. Our study suggests that the shift in spin chemical potentials generated by spin 
' accumulation produces an effective charge offsets. 

To demonstrate this charge offset effect, we apply a modified master equation calculation which takes into account 
the spin dependent charge states of the island. In this framework, the states are described by two parameters, Qj 
and Qi, denoting excess upspin charge and excess downspin charge, respectively. Comparing with spin-independent 
case, each primitive charge state \Q) is built up by spin charge states of Q = Qf + Qj,. Generally speaking, these spin 
dependent charge states are in non-equilibrium with the presence of spin accumulation. However, under limit of short 
energy relaxation time, the occupation distribution of a particular spin would form an equilibrium Fermi distribution, 
i.e. f(e — A*-f(n)- The numbers of net spin N = (Qj — Q[)/e is related to the spin chemical potentials and the density 
of states of the island p(e) as N = J dep(e) (f(e — /Xf) — f(e — pi)) . 

Although the electrostatic energy of each charge state is spin independent, the tunneling rate is spin dependent 
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in two ways: first, the effective tunneling resistances for the majority and minority spin tunneling processes are 
multiplied respectively by 2/(1 — P) and 2/(1 + P); second, the spin chemical potential shift's presence also modifies 
transition rates by changing the numbers of possible tunneling processes. Consequently, the master equation for each 
spin charge state, together with certain spin flipping transitions, reads 



dt 



= J2 {rT(*,il*±i,i)p*±ij + ri(i,i|»,i±i)p iii± i} 



;=s,d 



J2 {r , T (i±i,j|i,j) + ri(i,j±i|* ) i)}p tf + f^J , in 



Z=S.D 



in which i = Qi/e and j = Qi/e denote the numbers of upspins and downspins respectively, and T s (i' , is the 

tunneling rate for spin direction s(=j, j) in junction I (S for source and D for drain) from states to and 
Pij is the probability that the island is in state \i,j). In this equation, only sequential tunneling process is considered. 
By introducing an energy-independent spin relaxation time r s , the spin flipping transitions can be explicitly written 
as: 

dpij \ _ 1 [ U/S -U/6 \ 1 /7 1 + exp(-/3£7) (0 , 

dt ) si ~ r s \ 1 - exp (-/3[/) K+lj - 1 + 1 - exp (/?[/) P< " lj ' +1 J t s S 1 - exp {-(3U) PlJ 1 ' 

The first and second terms describe respectively the increase of probability due to up-to-down and down-to-up flip 
processes, and the third term is the decrease arising from the opposite processes. Here we assume that the up-to-down 
spin flipping rate is proportional to /(e — — f{ £ ~ Mi)) f° r electron with energy e. U = /if — /zj. = (i — j)6 is 

the chemical potential difference of upspins and downspins, and 6 ~ \/p is the energy level spacing, a constant for 
normal metal. For positive U, up-to-down spin flip is favorable, while for negative U , down-to-up spin flip dominates. 
Under these conditions the probabilities of major spin states with large \U\ are greatly reduced while suppressing the 
spin accumulation. 

Eq_(ffl) can be solved under the stationary condition given by dpij/dt — as described in the spin independent 
case.EHI Through a particular distribution of p^ , one can obtain the amount of spin accumulation, quantified as average 
chemical potential difference, U = -(i — j)$Pij, and the spin current for spin s tunneling through junction I, 

T l = e E + hj\i,j)-T l s (i - Pij- (3) 

ij 

If there is no spin flipping processes, the spin is conserved and the spin current passing through the source and drain 
junctions is the same. If the spin flips too quickly so as to completely destroy the spin accumulation, then the ratios 
between the two spin currents, for source and drain junctions, will be the same as polarization of source and 

drain electrodes, respectively. 

To gain an understanding about this phenomena, here, we perform a simulation using device parameters similar to 
those in experiments et aB: Rs = Ru =400kf2, Cs = CD=300aF, C g =0.8aF, Pg = Pu=0A. Because the resistances 
are much higher than quantum resistance Rq , the contribution due to higher order tunneling processes is negligible and 
only sequential tunneling process is included. We consider IVb characteristics and current-gate voltage dependences 
(IVg) with both parallel and anti-parallel alignment of leads under the no spin- flipping condition at a temperature 
of ksT/Ec = 0.1. In the parallel configuration, no particular feature is found because the ratio between two spin 
currents is simply the polarization 0.4, and total current is the same as that of the spin independent case. When 
the leads are in anti-parallel alignment, the calculation provides much more interesting results. The total current is 
smaller than that of the parallel case, and the high bias differential resistance is increased by a factor of 1/(1 — P 2 ) 
and shows a generic F-7N/F TMR effect. The differential TMR as a function of bias voltage also exemplifies expected 
oscillatory behaviors E 

The IV g characteristics shown in Fig. [l]a exhibit particularly different behaviors than from the parallel case. The 
upspin and downspin currents are only the same at V g = and V g = e/2C g . A closer inspection reveals that the 
peaks of two IV g curves with opposite spins shift with increasing bias voltage. This effect can be explained when 
we consider two separated spin transport channels. When the leads are in anti-parallel alignment, the source and 
drain resistances for a particular spin channel may differ by several times. This results in a step-like structure, called 
Coulomb staircase, in the IVb characteristics, and a distorted saw-tooth-like IV g modulation. i^Che Coulomb staircase 
effect can explain the TMR oscillation and the asymmetric gate dependence of spin currents.™ However, for a more 
rigorous study, it is necessary to include the spin entanglement term in the Hamiltonian, 2EcQiQi/e 2 . In fact, our 
calculations suggest that the results of the two methods differ especially at high bias voltages where both and Qi 
are large. 



2 



From the view point of spin accumulation, the raised upspin (lowered downspin) chemical potential effectively gives 
rise to a positive (negative) charge offset. At low bias voltage regime (V& < 2Ec/e), when V g is gradually raised from 
zero to maximal value(=e/2C s ), the electrical current increases due to suppression of Coulomb blockade. The spin 
accumulation is, in turn, enhanced by the increased current, and consequently there is a rise in both upspin chemical 
potential and the effective charge offset. Since within < V g < e/2C g region, the charge offset is an ascending function 
of V g , and the upspin current increases more rapidly than that of the zero spin accumulation. On the other hand, 
the increment of downspin current is less effective, because the downspin chemical potential decreases as V g increases. 
Therefore, at low bias regime, the IV g for upspin and downspin are tilted respectively toward lower and higher V g 
directions, and form saw-tooth-like IV g patterns. At bias voltages far beyond threshold (V& 3> 2Ec/e), the current is 
not much affected by V g and the spin chemical potential is less sensitive to V g . The shift in spin chemical potential 
(relative to the no spin accumulation case) AjU-fm increases (decreases) with 14. At Vj, = 6E c /e, A/i-m) is about 
±Ec, corresponding to a charge offset of about ±e. Consequently, as shown in Fig. |l]a, the upspin and downspin IV g 
characteristics shift by one period in respect with each other and differ from IV g characteristics at low bias voltages 
by half period. The total current is shown in Fig. [l)b, allowing a comparison with the experiments. 

For further investigation of the effect of applied gate voltage on two spin currents, we define a quantity describing 
the polarization of the tunneling current: Pj = (If— /^/(Jf+JjJ. Fig. || shows bias and gate voltage dependence of Pj. 
Such dependence suggests the possibility of using a ferromagnetic SET as a gate-controlled current polarizer. Because 
of large Pj values, the optimum operating regime is at low bias voltage. At eV g /C g = ±0.15 and eVb/Ec — 0.7, 
the current polarizations reach a maximum value of ±0.33. One can also explore the temperature dependence of the 
polarization current. There are two ways that the effects of temperature can enter, both leading to the destruction 
of current polarization. One is thermal activated charge fluctuation and the other is decrease of spin flip time. The 
former is automatically included in the master equation calculation and its effect is shown in Fig [|(a). At T — 0, the 
value of Pi can be as large as the polarization of the lead itself, while at fc^T > 0.5Ec, the gate charge effect becomes 
negligible. To evaluate the effect of the spin flip process, we assume an energy independent spin flipping time and 
an energy level spacing 5 of l^ieV in the island, and perform the calculations using the same device parameters as 
above. Fig ||(b) shows the current polarization for drain junction operating at Vb = 0.7Ec/e as a function of gate 
voltage at ksT/Ec = 0.1 under several spin flipping times. Clearly, when the spin flip time is short as compared 
with the tunneling time r t = e/I of approximately 10ns, the spin accumulation diminishes and Pj = —0.4, which is 
simply the polarization of the drain electrode. However, since the chemical potential is proportional to the island's 
density of states, the required spin flip time would be shorter for nanometer-sized normal-metal islands in which the 
level spacing is of the order of 10~ 8 ~ 10~ 9 eV, which is much smaller than the assumed value. It has been proposed 
that the criteria for spin accumulation is related to the tunneling resistance R t as t s S/Tl > Rt/ Rqu Our calculation 
results agree with this prediction. 

The co-tunneling processes, which are thus far not included in our calculations, can also give induce effective spin 
flipping. In the spin independent case, co-tunneling is a second order process that preserves the charge state but also 
produces current. In the Coulomb blockade regime, where sequential tunneling is suppressed, the current is mainly 
due to co-tunneling. For spin co-tunneling, there are spin-conserved and spin non-conserved processes. The latter is 
that a spin enters the island and an opposite spin leaves. The forward and backward co-tunneling rates r co , r co for 
F/N/F SET can be written aM: 
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The AE is the energy difference between initial and final states. For spin conserved co-tunneling, AE — ±eV, whereas 
for up-to-down and down-to-up co-tunneling AE — ±eV — U and AE = ±eV + U, respectively ('+' for forward and '-' 
for backward). Eg and E^ are energy changes of the tunneling processes Q — > Q ± e for source and drain junctions. 
Rs,eS and I?D,eff are the effective tunneling resistances for the source and drain junctions. Note that for anti-parallel 
configuration, Ps,effPD,eff product in Eq. (Q) is PsPd product multiplied by 4/ (1 — Ps)(l + Pd), 4/(1 + Ps)(l — Pd), 
4/(1 + Ps)(l + Pd) and 4/(1 — Ps)(l — Pd) for up-to-up, down-to-down, up-to-down and down-to-up forward co- 
tunneling events, respectively. 7 ± are decay rates for the final charge states Q ± e of the two processes and are given 
by 
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To investigate the co-tunneling spin flipping, we use tunneling resistances of 40k£7, which is closer to Rq, and leave 
other parameters unchanged. In the Coulomb blockade regime, the co-tunneling spin flipping rate can be as large as 
10 8 Hz, which is comparable to the tunneling rate of I/e ~ 10 9 Hz. However, different from the spin relaxation as 
discussed above, the spin flipping induced by co-tunneling for \U\ <C eVj does not have preferred direction. That is, it 
does not smear out the spin accumulation but rather only enhances spin fluctuation. Fig. |] shows the spin fluctuation 

SN = \]W - U 2 /5 as a function of gate voltage at Vb — 0.75Ec/e with and without consideration of co-tunneling 
processes. 

In summary, we proposed theoretically a gate-controlled polarized current in ferromagnetic single electron transis- 
tors. Using modified master equation formalism, we calculate spin sequential tunneling rates and spin current when 
spin accumulation is present. The spin accumulation-induced chemical potential shift behaves like charge offsets, 
producing interesting effects to the IV g characteristics. When the gate voltage is tuned away from the maximal and 
minimal gate voltages, the current passing through the junctions is polarized. The thermal fluctuation and spin 
flipping processes are both shown to suppress the effects from charge offset. The co-tunneling event provides effective 
spin flipping processes, but it only enhances spin fluctuation. Taking advantage of the gate dependence of polarized 
current in a ferromagnetic SET, one can utilize this system as a tunable current polarizer. 

We thank the members of the SET Group of the Institute of Physics, Academia Sinica, Taiwan for fruitful discus- 
sions. This research was partly funded by the Nation Science Council No. 89-2I12-M-001-033. 
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FIG. 1. Current-gate voltage dependences for (a) the spin currents from eVb = 0-5Ec (bottom) to 6.0-Ec(top) and (b) total 
current from eVb = 0.5Ec (bottom) to 4.0£c(top) with a step 0.5Ec for an F/N/F SET, whose parameters are described in 
the text. In (a) the dotted and solid curves represent respectively up-spin and down-spin currents. Notice that the peaks and 
valleys for the two spin currents appear at different gate voltages. At eVj = 4.5i?c the two currents shift half period while near 
eVb — G.OEc, they are the same. In (b) IV g curves shift with bias voltages: at eVt = 6.0Ec, it shifts by 0.5e. (c) The scheme 
of the F/N/F SET considered. The arrows indicate the magnetizations of electrodes. 



FIG. 2. (a) Polarization of current dependences on bias voltage at eV g /C g = 0.15(solid curve) and -0.25(dotted curve). At 
eVb/Ec = 0.7 and eV g /C g = ±0.15, Pi reaches a maximum value of ±0.33 while at eVb/Ec ~ 4.0 and eV g /C g — ±0.25, Pi 
has another maximum, which is approximately 0.04 with an opposite direction, (b) Polarization as a function of gate voltage 
at eV b /E c =0.5 to 4.0 with a step 0.5. 
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FIG. 3. (a) Gate voltage dependences of polarization at eVb/Ec=0.7 for temperatures ranging between and 0.5Ec/ks 
with a step 0.05 under no spin flip assumption. At T = 0, Pi cannot be denned for C g V g /e < 0.33 and C g V g /e > 0.67 since the 
current is zero within that range. The effect is more pronounced at lower temperatures, (b) The same dependence for the drain 
junction at eVb/Ec=0.7 and ksT/Ec = 0.1 with several spin flipping times r s , from bottom to top: 10ns, 100ns, 1/is, lO^is, 
100/is (solid curves), and oo (dotted curve). For r s = 10ns, Pi is fairly close to the P value of the drain electrode, i.e. no spin 
accumulation. For t s = lOOps, spin accumulation is almost the same as the non-flipping case and Pi increases dramatically. 

FIG. 4. Spin fluctuation SN v.s. gate voltage V g in the Coulomb blockade regime, Vb = O.lmV. The solid and dotted curves 
represent respectively the result for sequential tunneling with and without co-tunneling processes. The fluctuation is increased 
by about 20%. The inset shows the average spin number iV avc = U/S for the two cases. 
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